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Abstract 

Strong convergence results on tamed Euler schemes, which approximate stochastic differential 
equations with superlinearly growing drift coefficients that are locally one-sided Lipschitz contin- 
uous, are presented in this article. The diffusion coefficients are assumed to be locally Lipschitz 
continuous and have at most linear growth. Furthermore, the classical rate of convergence, 
i.e. one-half, for such schemes is recovered when the local Lipschitz continuity assumptions are 
replaced by global and, in addition, it is assumed that the drift coefficients satisfy polynomial 
Lipschitz continuity. 
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condition. 
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1 Introduction 

It is a well-known result that a stochastic differential equation (SDE) with a superlinearly grow- 
ing drift coefficient has a unique solution if, its drift and diffusion coefficients satisfy a suitable 
monotone growth condition, usually called coercivity, and a linear growth condition respectively 
along with some suitable local Lipschitz continuity conditions. One could refer to Krylov [7] and 
the references therein for more details. In addition, the almost sure convergence and convergence 
in probability of the corresponding (explicit) Euler approximations were proved by Gyongy [3]. 
However, Hutzenthaler, Jentzen and Kloeden [5] proved recently that the absolute moments of (the 
aforementioned) Euler approximations at a finite time could diverge to infinity as it was implied 
in Higham, Mao and Stuart [4]. In other words, the essential property of uniform integrability 
may not hold for such sequences and, thus, one could not obtain results on strong (in an £^'-sense) 
approximations although such results exist in the cases of almost sure convergence and convergence 
in probability. Furthermore, the introduction of accelerated Monte Carlo schemes provides a strong 
incentive for the study of strong approximations of SDEs since the latter is required for the efficient 
implementation of the former. More information on this topic can be found in Gile's seminal paper 
|T], Giles and Szpruch [2] and the references therein. 

Recently Hutzenthaler, Jentzen and Kloeden [6J introduced the notion of tamed Euler schemes 
in which the drift term is modified so that it is uniformly bounded. With such an approach. 



* Corresponding author. Email: s.sabanis@ed.ac.uk 



1 



they are able to prove that the tamed Euler scheme converges strongly (with rate one-half) to 
the exact solution of the SDE if the drift coefficient is globally one-sided Lipschitz continuous 
and has a derivative which grows (at most) polynomially. In addition, they assume that the 
diffusion coefficient of the SDE satisfies a global Lipschitz condition and it grows at most linearly. 
Furthermore, they offer a detailed review of the use of implicit schemes and a comparison between 
tamed Euler and implicit Euler schemes. Their comparison demonstrates that the implementation 
of implicit schemes requires significantly more computational effort than the tamed version. 

One however observes in [6j that for the "taming" of the drift coefficient, the term is used 
when it is known from classical literature that the standard strong convergence rate is one-half. 
In other words, one expects that the use of should be sufficient in order to control the drift 

coefficient and achieve strong convergence of the numerical scheme. In this article, a generalisation 
of the results of Hutzenthaler, Jentzen and Kloeden [6] is presented by using a variant of their 
tamed Euler method while a simpler proof is provided. It is proved that, even when global Lipschitz 
continuity conditions are replaced by local conditions, the tamed Euler scheme converges in to 
the exact solution of the SDE. Moreover, as a consequence of the aforementioned generalisation, 
the classical rate of convergence is obtained under the same assumptions as in [6] . 

We conclude this section by introducing some basic notation. The norm of a vector x G M*^ 
and the Hilbert-Schmidt norm of a matrix A € M'^^™ are respectively denoted by |x| and \A\. The 
transpose of a matrix A € R'^^™ is denoted by and the scalar product of two vectors x, y € is 
denoted by xy. The integer part of a real number x is denoted by [x]. Moreover, = £^(0, J-", P) 
denotes the space of random variables X with a norm \\X\\p := (E[|X|p])^^'' < oo for p > 0. 
Finally, B{V) denotes the a-algebra of Borel sets of topological spaces V. 



2 Main Result 

Let ($7, {/t}{f>o},-7^, P) be a filtered probability space satisfying the usual conditions, i.e. the fil- 
tration is increasing, right continuous and complete. Let {W^(i)}{t>o} be an m-dimensional Wiener 
martingale. Furthermore, it is assumed that b{t,x) and a{t,x) are S(M-i-) (8) jB(]R'^)-measurable 
functions which take values in MJ^ and W^^"^ respectively. For a fixed T > 0, let us consider an 
SDE given by 

dX{t) = b{t, X{t))dt + a{t, X{t))dW{t), Vt e [0, T], (1) 

with initial value X{Q) which is an almost surely finite J-Q-measurable random variable. 
For every n > 1, and any t G [0,T], the following tamed Euler scheme is defined 

dXn{t) = bn{t, Xn{Kn{t)))dt + a{t, Xn{Kn{t)))dW (t) (2) 

with initial value -'^(O), which is an J'o-nieasurable, almost surely finite, random variable and 
Kn := [nt]/t. Moreover, it is assumed that 

bn{t,x):=— — I ■.M t,x), (3) 
1 + n °'\b[t, x)\ 

for any f € [0, T], x G M'^ and a = 1/2. One then observes that 

|6„(t,x)| <min(n",|6(t,x)|). (4) 

Moreover, for every n > 1, one deduces immediately that 6„(t, x) is a /B(M-i-) (S> jB(M'^)-measurable 
function which take values in MJ^. 
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Let LP denote the set of nonnegative p-ih integrable functions on [0, T], i.e. to say if / G 
then 



T 



\fit)\Pdt<^. 



We make the following assumptions. 



A-1. There exists a positive constant K such that, 

2xb{t,x) V \a{t,x)\^ < K{1 + \x\^) 

for any t G [0, T] and x G M'^. 

A-2. For every R> 0, there exists a positive constant La such that, for any t G [0,T], 

{x-y)ib{t,x)-b{t,y)) V \a{t,x)-a{t,y)\'^ < Lr\x - 

for all \y\ < R. 

A-3. For every R> and p > 0, there exists A''^ G L^, such that 

sup \b{t,x)\ < NR{t) 

\x\<R 

for any t G [0, T]. 

A-4. For every p>0, E[|X(0)|p] < oo. 

Remark 1. Note that due for each n > 1, the norm of bn is a bounded function oft and x 
and, due to A-1, the norm of a has at most linear growth. This fact guarantees the existence of 
a unique solution to ([2]). Moreover, it guarantees that for each n > 1, all moments exist, each of 
which is bounded above by some value that dependents on n, i.e. 

sup E[\Xn{tW]<N (5) 
o<i<r 

for any p > 0, where N := A^(n,p, T, E|X(0)|'') is a positive constant. 

Theorem 1. Suppose A-1 - A-4 hold, then the tamed Euler scheme ([2]) converges to the true 
solution of SDE ([1]) in -sense, i.e. 



hm E 

n— >oo 



sup \X{t)-Xn{t)\P 
0<t<T 







for all p > 0. 

A-5. There exists positive constants I and L such that, for any t G [0,T], 
{x-y){b{t,x)-b{t,y)) V \a{t,x)-ait,y)\'^ < L\x - yf 

and 

\b{t,x) - b{t,y)\ < L(l + |xf + \y\^)\x - y\ 

for all X, y £ W^. 

Corollary 1. Suppose A-1 and A-3 A-5 hold, then the tamed Euler scheme 
true solution of SDE ([1]) in -sense with order 1/2, i.e. 



converges to the 



E 



sup \X{t) - Xn{t)\^ 
0<t<T 



for all p > 0, where C is a constant independent of n. 
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3 Moment bounds 

Lemma 1. Consider the tamed Euler scheme given by equation If for some p >2, 



and A-1 hold, then 



and 



sup sup E[|X„(t)|P] < oo 

n>l 0<t<T 



sup n^nit) - Xn{Kn{tW] < C^-^''' 
0<t<T 



sup E 

o<i<r 



I X„ (t ) - X„ ( (t ) ) n 6„ (i , X„ ( At„ (t ) ) ) |P 

where C is a positive constant independent of n. 
Proof. One immediately writes 



E 



J K„{t) J H„{t) 



{Kn{r)))dW{r)\^ 



for every t G [0, T], and thus, due to Holder's inequality, 



E\Xnit) - XninnitW <2P'^\t - Kn{t)r'E / 1 6„ (r, X„ (K„(r))) j^dr 

J Kn{t) 

ft 



+ 2P-1EI / a{r,Xn{Kn{r)))dW{r)\P. 

J Kn(t) 



One then observes that. 



2P-i|t _ Kn{t)r'E / \bn{r,Xn{Kn{r)))\Pdr < 2P-in(-i)P 

J Kn{t) 



and since sup„>x supf<2- E[|X„(t)|P] < oo and A-1 hold for some p > 2, then 

1P/2 



E 



Kn{t) 



a{r,Xn{Kn{r)))dW{r] 











<CE 


[/ 






. J K 



I Kn(t) 



[l + \Xn{Kn{rW)dr 



Substituting ([9]) and ([TO]) in ([8]) yields Furthermore, ([7]) holds trivially, since 



E 



\Xnit) - XnMtW\bnit,XnMt)m 



< E 



\Xn{t) - XnMt)W 



n 



<C 



is true due to ([6]) . 

Lemma 2. Suppose that A-1 and A-4 hold, then for some C := C(T, K, E[|X(0)p]), 



sup sup E 

n>l 0<t<T 



l^n(t)|' 



< c. 



Proof. Let us define 



/„(r):=E / {Xn{s) - Xn{Kn{s))bn{s,Xn{Kn{s)))ds 
L Jo 



Then, one calculates 

In{T) =E 



'-J0 ^ J Kn{s) 



6„(r,X„(«;„(r)))dr + ( a{r, Xn{Kn{r)))dW{r))bn{s, Xn{Kn{s)))ds 

J Kn{s) 



E 

k=0 



E 



bn{s,Xnik/n))K( / 6„(r,X„(fe/n))dr+ / air, Xnik/n))dWir 



ds 



E / bn{s,Xn{Kn{s))) / bn{r , Xn{Kn{r)))drds 



and thus 



|/n(r)|<E / / \bn{r,XnMr)))\drds 

'0 J K.n{s) 



(12) 



Furthermore, Ito's formula gives 

\Xn{t)\^ =|X(0)|2+2 f Xn{s)bn{s,XnMs)))ds+ f \a{s, Xn{Kn{s)))fds 

Jo Jo 

+ 2 / Xn{s)a{s,Xn{Kn{s)))dW{s) 

Jo 

<|X(0)|2 + 2/ X„(^j„(s))6„(s,X„(At„(s)))(is+ / |a(s,X„(K„,(s)))|2ds 
Jo Jo 

+ 2 [ {Xn{s)-Xn{Kn{s))bn{s,Xn{Kn{s)))ds + 2 f Xn{s)a{s , Xn{Kn{s)))dW {s) (13) 

Jo Jo 
and thus, due to A-1, §1 and i^, for any t G [0, T], 



E|X„(t)|2 <C7(1 +E|X(0)|2 +E / \Xn{Knis))\^ds) 

Jo 

<C(1 +E|X(0)p + /" sup E[X„(n)pds) 

Jo 0<M<S 

sup E!X„(n)p < C(l + E|X(0)p+ / sup EiX„(n)p(is) < oo 

<u<t Jo 0<u<s 



which implies, 



where the positive general constant C is independent of n. One then observes that the application 
of Gronwall's lemma yields 

sup E|X„(n)p < C 



0<u<T 



where C := C{T, K, K[\X{0)\^]). 



□ 



Lemma 3. Suppose that A-1 and A-4 holds, then for some C := C{p, T, K, ¥.[\X{0)\p]), 



E 



sup \X{t)\P 

0<t<T 



V supE 

n>l 



sup \Xn{t)\f 
0<t<T 



< C 



(14) 



for every p > 0. 

Proof. It is well known in the literature that the result 



E 



sup \X{t)\P 

0<t<T 



< C 



holds for every p > 0. One could consult, for example, Krylov (1980) for more details. 



In order to prove the second part of ()14l) . an inductive argument is used below. First, one 
chooses p = 2 and observes that due to Lemma [2] that 

sup sup E|X„(t)p < C 

n>l 0<t<T 

holds for some positive constant C := C(T,i^,E[|X(0)p]) which is independent of n. Thus, ([7|) 
from Lemma [T] holds true for p = 2 and one could use pSh to obtain the following estimate for 
q = 2p, i.e. q = A, 

E[sup \Xn{s)\'^] <C7(1 + E[|X(0)|'']+ / E\Xn{Kn{s))\'^ds 
0<s<t Jo 



+ / E[|X„(s)-X„(K„(s))r/2|6„,(s,X„(At„(s)))|^/2]ds 



+ E[sup I / Xn{u)aiu,XniKniu)))dW{u)\l/^]) 
0<s<t Jo 

and the application of the Burkholder-Davis-Gundy (BDG) inequality yields 



(15) 



E[sup |X„(s)|«] <C(1 + E[|X(0)|«] + / E[ sup |X„(n)|«]ds 

0<s<t Jo 0<u<s 

+ E[ r \Xn{s)\Ms,Xn{^n{sm'ds\''/^]), 

Jo 

where C denotes again a general constant which is independent of n. Thus, the application of 
Young's inequality yields 

E[sup |X„(s)|'?] <C(l + E[|X(0)r]+ ['e[ sup |X„(n)|^]ds + -^E[sup |X„(s)|'?] 

0<s<t Jo 0<u<s ^(-^ 0<s<t 

+ f \a{s,XMsm^ds\'i/^]) 
^ Jo 

which, due to A-1 and Holder's inequality implies that 

E[ sup \Xn{u)\^\ <C7(1 + E[[X(0)['?] + / E[ sup |X„(u)j'?]ds) < oo 

0<u<t Jo 0<u<s 



and thus the apphcation of Gronwall's lemma yields that 



E[ sup \Xn{t)\''] < C 
0<t<T 



(16) 



holds for some positive constant C := C(T, ET, E|X(0)|^) which is independent of n. Thus, ([7]) from 
Lemma [1] holds true for p = 4 and one could use (jlSp to obtain the estimate (jl6p for q = 2p, i.e. 
q = 8. Repeating the same procedure (by induction) one obtains the desired result (|14p . □ 



4 Proof of Main Result 

For every R > and n > 1, let us consider the stopping times 

Tn := inf{t > : \X{t)\ > R}, p^R := inf{t > : |X„(t)| > R} and u^r := tr A pnR. (17) 
Lemma 4. Suppose that A-3 holds, then for any R > and p > 



lim E[ / \b{s A VnR,Xn{Kn{s A VnR))) " bn{s A VnR, Xn{Hn{s A ynR)W]dt = 0. 

Proof. One immediately observes that for p >2 

E[ / \b{s A 1/„/J,X„(k„(s a l^ni?))) - bn{s A M„iJ,X„(K„(s A UnR)W] 

Jo 

'■^ \h{s ^UnR,Xn{^^n{s) AVnR))?'^ 



(18) 



+ n-°|6(s A UnR,Xn{tin{s A l^„ij)))|)P- 



< oo 



(19) 



due to A-3. Thus the application of the dominated convergence theorem yields the desired result. 

□ 

Proof of theorem Q]. Let p>2 and consider 

Xn{s) ■■= X{s A l^nR) " Xnis A f„i{). 

One observes immediately that 



E 



sup \X{t)-Xn{t)\P 

0<t<T 



<E 



sup \X{t) - X„(t)|PI| or p„fl<T} 

0<t<T ~ 



+ E 



sup \Xn{s)\P 
0<t<T 



(20) 



Then, by the application of Young's inequality for q > p and rj > one obtains 
E 



sup \X{t)-Xn{t)\P "iT[i<T or p„R<T} 
0<t<T 



<!^E 



q 

+ 



sup \X{t)-Xn{t)\'^ 
0<t<T 



p 



qj]' 



{tr < r or PnR < T) 



p/{q~p) 



|e 


"|^(Tij)r" 


+ E 


>n(PnK)r" 


1 




RP 







9^' 



<p/{<l-p)Rp' 



(21) 



7 



Furthermore, one defines 

/3n{s) := b{s A UnR,X{s A l^nR)) - 6n(s A X„ (k„ (s) A 

and 



Q„(s) := a{s A VnR,X{s A fnij)) - <t(s a l/nij,X„(Kn(s) A fnij)) 



to obtain 



with 



\Xn{t)\' 



Xn(s)/3n.(s) + |an(s)P ds + 2 Xn{s)anis)dW{s 



(22) 



Xn{s)^n{s) ={{X{s) - Xn{Kn{s)){b{s,X{s)) -b{s,Xn{Kn{sm 

+ {XMs))-Xnis))ib{s,X{s))-bn{s,Xn{Kn{sm} 

which implies, due to A-2 and A-3, 

Xn(s)/3„(s) < Jn{s) :=[i2LR + l)|xn(s)|' + i2LR + l)\Xnis) - 

+ \b{s,Xn{Kn{s))) - 6„(s, X„ (/«„ (s))) |^ 

+ 2NR{t)\Xn{s) - X„(aj„(s))|} I[,<,^^], 

whereas 

\an{s)\^ <2f>-'{Lf\xn{sW + Lf\Xn{s)-Xn{Knm''}Ths<.^„nV 



(23) 



(24) 



(25) 



Furthermore, in view of the above estimate (j25p . one observes that the apphcation of Young's 
inequality yields 



\Xnis)anis)\P/^ <2P/^{i2Lf + l)|Xn(s)r + LflXnis) " X„(K„(.))r} 



(26) 



Thus, from (|22p one obtains by applying Holder and BDG inequalities that 

E[ sup \Xn{uW] < CE f \\Jn{s)Y''^ + \an{s)\P + \Xn{s)an{s)f^ 



0<u<t 







ds, 



where C := C{p, T) is a positive constant, which in view of (j23|) . ()25p and (j26p yields 



0<u<t 



E[SUP |Xn(u)n<Cij/ ^Xn{s)\'' + nXn{s^UnR)-Xn{Kn{s)^VnR)\' 







ds, 



+ C 



E 



|5(s,X„(k;„(s))) - 6„(s,X„(fi;n(s)))rds 



+ / A^^/'(S)E|X„(S A - Xn{Kn{s) A l^ni?)!^/'^^ 



(27) 



where Cr := Cr{p, T, Lr) and (the redefined) C := C(p, T) are positive consants. The application 
of Grownwall inequality results in 

lim E[ sup |x„(t)n = 



for every R> due to Lemmas [T] and HI Finally, given an e > 0, one can choose r] small enough so 

e 

3' 



^2«C < i 



R large enough so 
and n large enough so 

to obtain due to ([20]) and ^ that 

E 

and thus prove the desired result. 



' -2C<i 



pI{Q'p)rp 3 



n sup \Xn{t)f] < - 
0<t<T <J 



sup \X{t)-Xn{t)f 
0<t<T 



< € 



□ 



5 Rate of Convergence 

First one observes that if A-4 and A-5 hold, then 

\b{t,x)\ < \bit,x)-b{t,0)\ + \bit,0)\ < L{1 + \x\^)\x\ + Noit) < N{t)il + \x\^+^) (28) 

for any t € [0, T] and x € M*^, where N{t) G L*' for any p > 0. 
Proof of Corollary [7]. First one rewrites ()23p in the following way 

Xn(s)/3n(s) ={{X{S) - Xn{s)){b{s, X{s)) - b{s, X„(s))) 

+ {X{s) - Xnis)){b{s,Xn{s)) - X„ (k„ (s) )) ) 

+ (X(s) -X„(s))(6(s,X„(K„(s))) -6„(s,X„(k„(s))))}I[,<,,^^] (29) 
and adjust accordingly, due to A-5, the upper bound J„ from (i24|) 

Xn{s)f3n{s) < Jn{s) : = {(L + l)\Xn{s)\^ + L\1 + |X„(s)|' + |X„ (aC„ (s)) |' )' |X„(s) - X„(aC„(s))|2 



+ \b{s,XniKn{s))) -bn{s,XnMs)))\^}^s 



(30) 



and, thus, the last term of (j27p is replaced by 



£{t) :=E 



C{1 + |X„(S)|'P + \Xn{Kn{s))\'^)\Xn{s) - X^^sW 



ds 



which is estimated from above by 



£{t) <C I (VE|X„(s a VnR) - Xn{Hn{s) A l^„ij)j2p)ds 
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due to Holder's inequality and (jl4p . Note that the general constant C is independent of t and n. 
In view of Lemma [H one deduces that 

sup Sit) < Cn-P/\ (31) 

0<t<T 



Furthermore, ([T9]), ([11]) and ([28]) imply that 
rT 

E[ / \b{s A Z^„K,^n(«:n(s A UnR))) - 6„(s A UnR,Xn(,Kn{s A J^„K)))r] < Cn^^/^ (32) 







which along with ([6]), (i3T]) and (1321) result in 



E[ sup |xn(t)n <C^^'^/' (33) 

0<t<T 



due to ()27p . Finally, one chooses rj = n 2 ^ R = n'^ii-p) , q > p > 2, to obtain the desired result due 
to ([20]), ^ and (jajl . 

□ 
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